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Abstract 

The problem of community detection is relevant in many disciplines of science and modularity 
optimization is the widely accepted method for this purpose. It has recently been shown that 
this approach presents a resolution limit by which it is not possible to detect communities with 
sizes smaller than a threshold which depends on the network size. Moreover, it might happen that 
the communities resulting from such an approach do not satisfy the usual qualitative definition of 
commune, i.e., nodes in a commune are more connected among themselves than to nodes outside the 
commune. In this article we introduce a new method for community detection in complex networks. 
We define new merit factors based on the weak and strong community definitions formulated by 
Radicchi et al (Proc. Nat. Acad. Sci. USA 101, 2658-2663 (2004)) and we show that this local 
definitions avoid the resolution limit problem found in the modularity optimization approach. 
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I. INTRODUCTION 



The problem of community detection in complex networks has recently attracted the 
attention of researchers in different areas of scientific knowledge. This is due to the fact 
that it is a common practice to represent some complex systems as networks constituted by 
interconnected nodes. 

A network G is defined by a set of nodes {n} (rii, n 2 ,...,n n .), and a set of links {1} 
(li2,lu---,hm)- A link ly denotes a relation between node n« and node rij. Depending on 
the possible values of Uj the resulting network can be of two types. If Uj can only have the 
values 1 or 0, we will call the network unweighted, on the other hand, a network will be 
defined as weighted if can attain values different from or 1, thus indicating that the 
relation between nodes is also characterized by a given strength. In this work we will focus 
on unweighted networks. We will assume that for every node 7ij there exists at least another 
node rij such that kj is different from 0, moreover we will consider networks such that for 
every conceivable pair of nodes there will be a path (i.e. a sequence of links {hjljkhm---}) 
joining them, in such a case we say that we are dealing with connected networks. We will 
consider that the links are undirected i.e. kj = Iji. Further on, we will focus on sparse 
networks for which the number of links in {1}, L, is much less than the maximum possible 
number of links, L max given by L max = N(N — l)/2 , with N the total number of nodes in 

W • 

Generally, complex networks contain a large number of nodes and links, and it is often 
possible to decompose them into subgraphs called communities, selected according to a given 
criterion. A community is usually defined, qualitatively, as a subgraph of the network whose 
nodes are more connected among them than to nodes outside the subgraph 

Community detection has a wide range of applications. The partition of a network in 
communities might allow us to find an specific function naturally assigned to each commu- 
nity, as for example in the case of metabolic networks [3(]. On the other hand, community 
detection can help us to identify social groups in a social network or can be used to perform 
a coarse-graining reduction of the network to simplify subsequent analysis {4]. 

There are many methods to decompose a network into communities, but the widely 
adopted in recent years is the one proposed by Newman and Girvan |2|. These authors 
define a merit factor named modularity (Qn) that quantifies the quality of a given m- 
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subgraphs partition {Cj}i<j< m of the graph G, where \JT =1 Cj = G with Cif]Cj = if 
i ^ j. This quantity measures the difference between the actual fraction of internal links 
in each subgraph with respect to the expected value of the same quantity, if nodes in the 
network are randomly connected keeping the degree of each one fixed. The best partition of 
the network is taken as the one that maximizes the modularity Qn, in this way, the network 
partition problem is turned into an optimization one. 

Modularity optimization is a hard problem due to the fact that the number of possible 
partitions of a network increases at least exponentially with its size. Indeed, it has recently 
been proven that this problem is NP-complete |5| , and then, there is not a correct polynomial- 
time algorithm to solve it for networks of any size. Many optimization algorithms have been 
developed, like simulated annealing si, [f3], extrema optimization 7] and spectral division Q], 
but all of them can only give an approximation to the optimum partition for large networks. 

In this work, we do not introduce a new Qn optimization algorithm, but we propose new 
merit factors for the calculation of the partition of networks into communities. 

The modularity Qn is a non-local community definition in the sense that it is necessary to 
know general characteristics of the whole network in order to decide if a given subgraph of the 
network is a community. In a recent paper [lj Fortunato & Barthelemy have shown that this 
non-local character imposes a resolution limit, by which the minimal community size that 
can be detected, by modularity optimization, depends on global network parameters. Then, 
Qn optimization is not able to detect communities of size smaller than a given threshold. 

In this work we use the weak and strong community definitions proposed by Radicchi et 
al in 1(3] . We emphasize their local character and introduce new merit factors to evaluate 
the quality of a given partition of a network based on these quantitative definitions of 
community. Then, we implement an optimization method in the spirit of simulated annealing 
, in order to analyze different networks using these new merit factors and to show the 
characteristics of our approach. Finally, we show that the resolution limit problem does not 
appear in our approach. 

The paper is organized as follows. In section [III we review the definitions introduced 
in [lo| and compare them with the modularity Qn, we also analyze the meaning of this 
last quantity. In section III Al we define the community strength S in the strong and weak 
sense and we introduce the associated merit factors. In section imi we apply our method to 
different, well known, networks. In section [IV] we analyze the resolution limit problem for 
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our approach. Finally, conclusions are drawn in section |VJ 



II. COMMUNITY DEFINITIONS 

When thinking about communities in networks we have in mind a qualitative community 
definition: a community is a group of nodes in which the number of internal links, connecting 
nodes within the group, is larger than the number of external ones. In order to formalize 
this qualitative criterion we consider a graph G containing N nodes, with ki the degree of 
node i £ G. If C is a subgraph of G with kf 1 and k° ut the number of links of node % G C 
that connect it to nodes inside and outside of C respectively. There are two quantitative 



community definitions introduced by Radicchi et al [1 Ql j : 



i) Community in strong sense: C is a community in the strong sense if: 

k m > k out y ieC ^ 

ii) Community in weak sense: C is a community in weak sense if: 

K > E k T l (2) 



In words: a subgraph C C G will be a community in the strong sense if each of its nodes 
has more links connecting it with nodes in C than those that connect it with other nodes 
not belonging to C. In a similar way, C C G will be a community in the weak sense if the 
sum of the number of links that interconnect nodes inside C is larger than the sum of all 
links that connect nodes in C with nodes not belonging to C. These community definitions 
are simple, intuitive and local: given a subgraph CcGwe can decide if it constitutes a 
community, in either strong or weak sense, without knowledge of the entire structure of G. 

In order to compare the previous approach with the one proposed in j^, we briefly review 
the definition and meaning of the modularity Qn- Given a m-subgraphs partition {Cj}i<j< m 
of the graph G, where Uj=i Q = G> the mathematical expression of Qn is : 

m |" ; / j \ 2 



h ( dj 
L \2L 



(3) 



where U denotes the total number of internal links for subgraph Gi C G and &{ = Xljec %> 
and L — | Ylij^c k j is the total number of links in G. 



The term U/L in Eq. [3] denotes the actual fraction of internal links in subgraph Cj, while 
di/2L can be interpreted as the probability of a link to be connected to a node in subgraph 
Cj. Then, (di/2L) 2 is the expected fraction of links within subgraph Cj when all nodes in 
G are randomly connected keeping the degree of the nodes fixed. This last ideal random 
picture is used to compare with the actual one because it is assumed that corresponds to a 
situation with no communities (although it was shown in ll| that random networks may 
have a community structure). 

As already mentioned, the modularity Qn was conceived as a measure of the goodness 
of a given partition of the network. Then, the bigger Qn is, the better the partition is. 
We should notice that this merit factor implies, in turn, a community definition (which 
does not necessarily corresponds to the intuitive one stated above): a subgraph Cj will be 
a community if the actual number of links that connects nodes in Cj is bigger than the 
expected one when all nodes in the network are randomly connected, this is to say, when 
li/L — (di/2L) 2 > 0. Clearly this last condition depends on the global parameter L , then, 
we say that the community definition associated with Qn is non-local. 

In what follows, we will introduce new merit factors associated with the weak and strong 
community definitions. 



A. Merit factors for weak and strong 
community definitions. Community strength 

Given a graph G and a m-subgraphs partition {Cj}i<j< m , where each subgraph Cj C G 
constitutes a community according to any of the local definitions mentioned in the previous 
section, we want to define a quantity that measures the "quality" of each of the resulting 
communities. In the context of the above mentioned local framework, this quantity must 
only depend on the local characteristics of the subgraph Cj. Therefore, our analysis must be 
circumscribed to nodes and links belonging to Cj and external links that connect nodes in 
Cj to nodes outside Cj. Following the weak and strong definitions of community, the more 
internal links a community has, with respect to the external ones, the "stronger" it will be. 
If ki = k] n + k° ut is the degree of node % G Cj, where k\ n and k° ut are the number of internal 
and external links for node i, we define the "community strength" (S), that measures the 
normalized difference between internal and external links for nodes in Cj\ 
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h,in uout 

ieCj 

were L(Cj) = \ Y^iec Then, — 1 < S(Cj) < 1, and it achieves its maximum value 1 
when kf"* = Vi G Cj. 

The definition of S(Cj) according to Eq. H]is valid for unweighted networks. In the case 
of weighted links, we have to interpret fcj as the sum of the weights of the links that connect 
to node i, for both k\ n and kf ut . 

Now we introduce a merit factor Qw for the weak community definition as the sum of 
S(Cj) over all subgraphs Cj C G: 



m m uin uout 

3=1 j=l ieCj K 3 > 



with the constraint that each subgraph Cj C {Cj}i<j<m must satisfy the weak community 
definition i.e. 



S(C 5 ) > VC, C {C,} 1 <,< m (6) 

As in the case of Qn- the bigger Q w is, the better the m-subgraphs partition {Cj}i<j< m 
of G will be, in the sense of weak community definition. Then, it is possible to implement 
the optimization algorithms developed for for this new merit factor Qw, 

In Eq. El Qw — 1 when all the network constitutes a single community. If as a result of 
the optimization process the maximum value obtained for Qw is precisely Qw = 1 and we 
get a single community, then the best partition of the graph corresponds to no partition. 
However, it is possible that one could get S(Cj) > WCj C {Cj}i<j<m for a given m- 
subgraphs partition, with m > 1 but with < Qw < 1- The resulting community structure 
of network would correspond to a suboptimal partition. 

In the same spirit we now define a merit factor Q$ according to the strong community 
definition: 

m m h,in uout 

3=1 j i r, 1 i> 

with the constraint 

(Af - k° ut ) > Vz G Cj (8) 
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Now, our definition of optimal partition can be stated in the following way: 
Definition: the optimal m-subgraphs partition {Cj}i<j<m of a graph G in the Strong 

(Weak) sense is that one with maximal merit factor Qg ( Qw )■ 

In next section we will show some examples of the application of this new merit factors 

in network partition problems. 



III. EXAMPLES 



In all examples presented in this section we have used an optimization algorithm based 
on simulated annealing, described in [6|, but for our new merit factors. The optimization 
can be performed in two ways. In the first one the total number of communes is left as a 
free parameter and as a consequence the final number of communes is determined by the 
simulated annealing process. In the second one, the number of communes is taken as an 
extra constraint. We will always use the first approach unless it is explicitly stated that 
the number of communes is fixed. This last methodology might be used when the optimal 
number of communes is already known from the experiments as in the case of the Zachary 
network. 



A. Zachary's karate club network. 



We will begin with a typical case: Zachary's Karate Club [12|, that has turned into an 
unavoidable example in publications about community structure. This network represents 
the relationships between members of a karate club at a University in the 1970s, and it 
has been shown that it has a strong community structure in many previous studies 
Applying the optimization algorithm for the weak community definition merit factor Qw-, we 
obtained, for the unweighted version of Zachary network, a partition into three communities 
of sizes: 17 (Ci), 12 (C 2 ) and 5 (C 3 ) nodes, with a value of Q w = 1.792 (Fig. Q}. When the 
number of communities was constrained to two, we obtained two communities of 17 nodes 
each one, with Qw = 1.487 (circles and squares in Fig. CD). This partition corresponds to 
the one observed by Zachary with the exception of node 10 that appear misclassified. 

With this analysis we can know, in addition, the strength S(Cj) of each community Cj 
in the network. For the best partition of the Zachary network into three communities of 17 
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FIG. 1: (Color online) Best partition for Zachary network. Squares and circles denote the two 
communities obtained with our approach when the number of communities is fixed to two. This 
partition correspond to the one consigned by Zachary in [2], with exception of node 10 that appear 
misclassified. 

(d), 12 (C a ) and 5 (C 3 ) nodes, we have: S(d) = 0.744 , S(C 2 ) = 0.548 and S(C 3 ) = 0.5, 
with C\ as the strongest community. On the other hand, the partition into two communities, 
is composed by two strong communities of 17 nodes each one, with S(Cj) = 0.744 for each 
one. 

When we perform the community analysis using the strong community merit factor Q$, 
we obtained two communities: C\ with 29 nodes [S{C\) = 0.943) and C 2 with 5 nodes 
(iS^Cy = 0.5). In Fig. Q] can be observed that node 10 has one internal and one external 
link and this situation is not be allowed in the strong community definition. For this reason, 
the communities with 17 and 12 nodes are joined together. 

B. Star network. 

Another testing example is the star network of Fig. [2], consisting of two interconnected 
stars with 11 nodes each one. In the weak community picture we obtain for the optimal 
community structure a partition into two communities of 11 nodes each with a value of 
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FIG. 2: (Color online) Community structure for star network obtained in our approach. The two 
communities are distinguished with squares and circles. For Qn optimization we find an additional 
community containing nodes 4 and 21 with strength S = 0. 

Qw = 1.652. The corresponding strengths are S(C\) = 0.833 (circles in Fig. [2]) and 
S(C2) = 0.818 (squares in Fig. [2]). The difference in strength between the two communities 
is ascribed to the extra link that connects nodes 2 and 3. 

We must notice that, in the same context, there is another partition with the same value 
of Qw but composed by two communities with 12 and 10 nodes. This happens when node 
4 is moved from one community to the other in Fig. [2j 

No partition was obtained when we use the strong community merit factor Qs- This is 
due the fact that nodes 4 and 21 are singly connected and then the condition kf 1 — k° ut > 
is not satisfied. 

When Qn optimization was implemented, we obtained three communities: C\ and C*2 
with 10 nodes each and C3 with 2 nodes, including nodes 4 and 21, with strength S{C^) = 
which does not satisfy any of the quantitative community definitions reviewed in [Til 
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FIG. 3: (Color online) Ring network. This symmetric network does not present an intrinsic 
community structure. However, almost all approaches for community detection find a community 
structure for this example, but the identities of the nodes within each community change when the 
detection process is repeated. 

C. Ring network. 

Another example is the ring network of Fig. [3] with 20 nodes and k = 6. This network 
can not have a community structure due to its symmetry. However, we have obtained two 
communities of 10 nodes each one, with strength S = 0.6 by means of the weak merit 
factor optimization. We must notice here, that the found communities are not unique, that 
is to say, on having applied repeatedly the algorithm different communities of the same 
size but involving different sequences of node indexes are obtained. This is an evidence 
of the absence of an underlying community structure. This unsatisfactory result is also 
obtained when optimizing the Qn merit factor, but in this case the optimal partition is into 
3 communities, two of them with 7 nodes each, and the third one with 6 nodes 

On the other hand, when we run the optimization algorithm with the strong community 
condition, no partition is obtained. This is true for all ring networks because an hypothetical 
frontier node, with kf 1 = k° ut , will not satisfy the strong community condition. 
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D. The bottlenose dolphins network 



Another social network which has attracted considerable interest is the one corresponding 
to the bottlenose dolphins network, which has been fully analyzed in [14] (see also {2]). This 
small social network is composed by 62 nodes and it is known to consist of two communities 
of sizes 41 and 21 nodes each. Following the approach proposed in this work we first analyzes 
this network applying the Qn analysis in our simulated annealing approach. The result of 
this analysis is the partition of the network into four communities composed by 21, 16, 
13, and 12 nodes each. When we performed the optimization of the Weak community 
definition we obtained five communities of 20, 12, 11, 10 and 9 nodes each. Finally when 
the dolphin network is analyzed in terms of the strong community definition we obtained 
the actual partition, as observed experimentally, in two communities of 41 and 21 nodes 
each. These last two results are displayed in Fig. HI In this figure we show the two 
communities according to the Strong community definition as circles (41 nodes community) 
and as squares (21 nodes community). The corresponding analysis according to the Weak 
community definition further divides the previous two communities and are denoted by the 
different shades of gray (see caption for details) in the figure. It should be noted at this point 
that when the optimization of Qw is performed with the extra constraint that the number 
of communes is 2 we obtain the same community structure as observed experimentally. 



E. Computational generated test network. 

We conclude our short list of examples with the analysis of computer generated graphs 
which have a community structure. A word of caution should be raised at this point, because 
the communities built into these graphs are usually of uncertain nature and being purely 
theoretical, the assumed community structure cannot be "verified experimentally" as in the 
case of the Zachary Karate Club or the bottlenose dolphins case. 



1. Methodology I 

In this case we use the method proposed and analyzed in 0, Q. We take a 128 nodes 
graph G divided in four modules Cj=i t ,,,^ of 32 nodes each and with nodes degree ki = 
k] n + k° ut = 16 Vi G G. We have earlier defined k™ (k° ut ) as the number of links that 
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FIG. 4: (Color online) Bottlenose dolphin network. This network has a size of 62 nodes and it 
is known from direct observation that it has two communities. In this figure squares and circles 
denote the communities detected by our Strong community approach and the colors (shades of gray 
or colors on line) show the results of the Weak community approach. Notice that the optimization 
according to Qw merely subdivides the communities obtained through Qs optimization. 

connect node i G Cj to another node in (out of) Cj. When k° ut is varied between to 16, 
G go from a strong communality graph to a quasi-random one. 

In the framework of weak and strong merit factors optimization the expected four com- 
munities partition was obtained for < k° ut < 7, where Qw < 1 only for k° ut = 7. When 
k° ut > 8 the mean number of external links is bigger or equal to internal ones in each module 
and no partition was obtained. 
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2. Methodology II 



We now use the formalism introduced in 



15] , in which an algorithm for generating a class 



of benchmark graphs, that account for the heterogeneity in the distributions of node degrees 
and of community sizes was devised. The aim of this algorithm is to built a graph with more 
or less well defined community structure (in the caption of Figure 5 in this paper a reference 
is made to "communities in the strong sense" which is not the case for this algorithm as is 
easily verified). It is assumed that both the degree and the community size distributions are 
power laws, with exponents 7 and (3, respectively. The number of nodes is N, the average 
degree is < k >. One more parameter characterizing this model is the mixing parameter fi. 
Each node shares a fraction 1 — fi of its links with the other nodes of its community and a 
fraction /x with the other nodes of the network. 

We have generated graphs according to this algorithm and we have analyzed them using 
the Girvan-Newman definition of community using our global optimization approach [6] and 
the one proposed in this work. The parameters defining the graphs were chosen to be: 
7 = 2.5, P = 1.5, 0.1 < fi < 0.6. The size of the graphs was fixed in 300 nodes, and the 
mean degree in 8 (maximum degree=30). 

The results of such a calculation are summarized in Fig. [5] (see captions for details). 
The quality of the partitions obtained with the recognition algorithms with respect to the 
communities established by the algorithm of Lancichinetti et al. is measured by means of 
the so called Normalized Mutual Information [16]. According to the results displayed in Fig. 
our merit factors definitions outperforms the Girvan-Newman approach. 



IV. RESOLUTION LIMIT PROBLEM 

In a recent paper |l[ S. Fortunato and M. Barthelemy showed that modularity Qn op- 
timization fails to detect communities smaller than a certain threshold which depends on 
global parameters of the network under study, as is the case of the total number of links in 
the network. 

Following [l| we can define a community from Eq. [3j in the framework of modularity 
Qn, as a subgraph Cj C G that satisfy: 
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0,1 0,2 0,3 0,4 0,5 

Mixing parameter 



FIG. 5: (Color on line) Computer generated test network. We have generated graphs according 
to the formalism proposed by Lancichinetti et al. The results of our calculation are displayed by 
the full circles (in red) while the results according to the Girvan-Newman approach are denoted 
by square (in blue). 



i-Qi) 2 > o (9) 

This expression can be interpreted as the community strength in this framework. We can 
write di = 2/j + l° ut , where U and l° ut denote the number of internal and external links for 
subgraph d, and write l° ut as a fraction of internal links l° ut = all with a > 0. Then, from 
Eq. [9j the following condition for the community size U is obtained: 

<< < uTW (10) 



y shows that the community 
l| the authors show that this 



The dependency on total number of links L in Eq. [9] clear 
definition, in the context of modularity Qn, is non-local. In 
non-locality is the origin of limit resolution problem. 

On the other hand, in our approach the community definition is strictly local. Then, we 
can decide if a subgraph is a community without regards to the size of the entire network. 
To illustrate this conclusion we work out an example introduced in 



14 



FIG. 6: Ring of cliques. Each circle represents a clique which is a totally connected subgraph with 
n nodes. 



Let us suppose a ring of totally connected subgraphs (from now on cliques) in Fig. El 
Each subgraph has n nodes connected by n(n — l)/2 internal links and two external ones, 
and we have M of this subgraphs with a total number of links L = M(n(n — l)/2 + 1). The 
optimal partition for strong and weak community definition frameworks is the natural one: 
each clique constitutes a single community. This can be easily shown by analyzing another 
alternative partition in which each community Cj contains W > 2 cliques. 

First, we calculate the strength S(Cj) for one of this subgraphs containing W cliques: 

b[Cj) ~ W(n-l)n + 2W W 

Now, we want to compare the result of Eq. [11] with the total strength of the same subgraph 
Cj when each clique is taken as a single community. In this case, the strength of Cj is given 
by W times the strength of one clique: 

v 3J n(n- 1) + 2 K J 

Then, it is straightforward to see that S*(Cj) > S(Cj) is equivalent to: 

Win - 2)[{W - l)(n - 1) + 2W] - 2{nW - 2) > (13) 

We have said that each clique constitutes a single community in the strong (and then, 
also in the weak) community definition, therefore n > 3. Then, the condition of Eq. [13] is 
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satisfied for all W > 2. This is to say that the optimal partition in our approach is that 
one for which each clique constitutes a single community. This is not the general case in 
the modularity Qn framework. It was showed in [1] that, due to resolution limit problem, 
partitions in communities with two or more cliques can give larger values of Qn than with 
single clique communities. 



V. CONCLUSIONS 



In this work we have proposed new merit factors to recognize communities in networks. 
These merit factors are more realistic than the ones currently in use in the literature because 
they strictly adhere to what a community is expected to be, i.e., a subset of nodes which 
are more connected among themselves than to the rest of the network under consideration. 

We started by putting forward this qualitative definition of a community and then we 
reviewed the meaning of the quite popular measure of the quality of a given partition known 
as the modularity Qn- As we have discussed above, the community definition associated 
to this quantity is non-local and does not necessarily corresponds to the aforementioned 
qualitative definition. One of the consequences of the non-local character intrinsic to this 
quantity is the limit resolution problem as stated in jj. 

In order to recognize communities in networks that strictly adhere to the qualitative 
definition, we have used (following [loj]) two local community definitions: the weak one and 
the strong one . In order to use this definitions to recognize communities we have developed 
a criteria to quantify the strength of a community (S). Afterwards, we have defined two 
merit factors associated with S which we named Qs and Qw- As with Qn the problem 
of recognizing communities in a network is mapped onto an optimization problem, i.e., the 
communities in a network are the elements of the partition which maximizes Qs or Qw 
We have performed the optimization of these merit factors on some standard networks by 
implementing an algorithm in the spirit of simulated annealing. The limit resolution intrinsic 
to the Qn definition is not present in our approach. 

It is worth noticing at this point that the solution to the detection of communities in the 
strong sense is also a solution in the weak sense but not necessarily optimal. On the other 
hand, the converse is generally not true as we stated in section IIHI 

The strong community definition tends to give larger communities because of its inability 
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to deal with nodes that are equally shared by two highly connected subgraphs, but on 
the other hand has the nice property that it is the only one that gives no partition for 
symmetric string networks and also solves the problem of the bottle nose dolphins network 
exactly, without constraints in the number of communities. 
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